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CONTRACT  STRESS  ANALYSIS 
OF  CERAMIC-TO-METAL  INTERFACES 
FINAL  REPORT 


1.0  INTRODUCTION 

This  Pinal  Report  presents  the  technical  progress  completed  for 

I 

Contract  N00014-78-C-0547,  "Contact  Stress  Analysis  of  Ceramic-to- 
Metal  Interfaces." 

The  Navy  Scientific  Officer  is  Mr.  Bob  Pohanka. 

The  objective  of  this  proposed  8-montb  program  was  to  conduct 
analysis,  specimen  testing,  and  data  correlation  to  provide  an 
improved  understanding  of  the  local  contact  conditions  that  prevail  at 
an  interface  between  ceramic  and  raetal  components  for  gas  turbine 
engines.  The  program  was  specifically  directed  to  study  contact 
stresses  at  the  interface  between  in&er:ed  ceramic  turbine  blades 
(hot-pressed  silicon  nitride)  an^'  a  metal  rotor,  but  the  method  of 
analysis,  results,  and  conclus  ona  -Iso  provided  a  better  under¬ 
standing  of  contact  stresses  at  ceramic-to-metal  and  ceramic-to- 
ceramic  interfaces  for  static  components. 

A  finite-element  stress  analysis  procedure  and  structural  evalua¬ 
tion  technique  consistent  with  the  statistical  nature  of  ceramic  mate¬ 
rials  was  generated.  Specimen  testing  was  conducted  to  obtain  vali¬ 
dation  of  these  techniques  ''d  provide  insight  as  to  possible  modifi¬ 
cations  of  stress  simulation  or  fracture  prediction  criteria.  The 
objective  was  to  develc^^  a  design  methodology,  which  will  Improve  on 
current  design  methods,  thus  permitting  the  design  of  interface  con¬ 
figurations  in  which  strengths  more  indicative  of  the  inherent 
strength  of  ceramic  materials  can  be  realized.,. 

N. 

This  four-task  effort  was  subdivided  as  foIlowsTv 
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Task  I 
Task  II 
Task  III 
Task  IV 


Boundary  Condition  Study 

Finite-Element  Model  Development  and  Analysis 
Specimen  Testing  and  Model  Validation 
Reporting 
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2.0  TECHNICAL  PROGRESS  SUMMARY 

Task  I  of  the  program,  Boundary  Condition  Study,  has  been  com¬ 
pleted.  A  simplified  model,  a  cylinder  on  a  plane,  was  devised  which 
could  be  analyzed  with  both  a  closed-form  solution  and  a  finite- 
element  computer  program.  This  model  is  also  amenable  to  specimen 
testing. 

Task  II  of  the  program,  Finite-Element  Model  Development  and 
Analysis,  has  been  completed.  Significant  progress  has  been  made  in 
modeling  the  contact  zone  of  the  ceramic-to-metal  interface  region.  A 
cylinder  in  contact  with  a  plane  was  chosen  for  the  analytical 
approach  because  it  not  only  approximates  the  situation  of  ceramic 
blade  to  rotor  attachments,  but  it  provides  closed-form  and  finite- 
element  solutions  for  both  the  non-friction  and  friction  cases. 

Task  III,  Specimen  Testing  and  Model  Validation,  has  been  com¬ 
pleted.  A  test  fixture  for  the  verification  testing  has  been  designed 
and  fabricated  to  apply  a  tensile  load  to  the  ceramic  specimens  while 
applying  a  compressive  normal  contact  load.  Various  ratios  between 
compressive  and  tensile  loads  have  been  investigated  with  and  without 
compliant  layers.  Both  HS25  and  platinum  compliant  material  have  been 
investigated. 

Task  IV,  Reporting,  has  been  completed  with  the  submission  of  the 
Final  Report,  which  documents  the  work  completed  within  this  8-ffiOnth 
program. 

2.1  Analytical  Approach  and  Accomplishments 

2ooa  finite-element  modeling  techniques  have  been  used  success¬ 
fully  to  calculate  stresses  in  the  interface  between  bodies  in  con¬ 
tact.  Good  agreement  has  been  achieved  between  closed-form  solutions, 
References  1  and  2,  and  finite-element  models  for  a  cylinder  contact¬ 
ing  a  semi-infinite  solid  with  and  without  friction  forces.  These 
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particular  solutions  approximate  a  ceramic  blade-to~metallic  disk 
attachment. 

As  previously  reported,  a  finite-element  approach  was  used  to 
analyze  the  frictionless  case  for  a  ceramic  cylinder  and  a  ceramic 

plate.  The  analysis  employed  zoom  modeling  to  provide  a  sufficient 

number  of  nodes  in  the  region  of  contact.  The  original  model, 
Figure  1,  provided  boundary  conditions  for  a  more  refined  model  of  the 
contact  zone,  Figure  2.  This  model  did  not  have  sufficient  refinement 
to  give  adequate  representation  of  the  contact  zone.  Consequently, 
results  for  this  model  provided  input  for  a  more  refined  finite- 

element  mesh.  As  before,  this  refined  model  still  did  not  have  suf¬ 

ficient  nodes  in  the  interface  region  to  completely  depict  the  stress 
distribution,  so  a  final  zoom  model  was  generated  which  had  thirteen 
nodes  in  the  contact  zone.  Figure  3.  The  relative  sizes  of  the  four 
finite-element  models  ate  shown  in  Figure  4.  Figure  5  compares  the 
stresses  in  the  direction  normal  to  the  plane  for  the  four  finite- 
element  models  and  a  classical  closed-form  solution  from  Reference  1, 
Good  agreement  was  also  achieved  for  maximum  shear  stresses  and  trans¬ 
verse  stresses. 

Contact  between  dissimilar  materials  was  also  studicnS.  An  analy¬ 
sis  similar  to  that  mentioned  above  was  completed  for  an  aluminum 

cylinder,  E  «  7.10  x  10  N/cm  ,  and  a  ceramic  plate, 
7  2 

E  “  20,96  X  10  N/cm  ,  Results  from  the  zoom  finite-element  solution 
agreed  well  with  the  closed-form  results,  Figure  6. 

The  transverse  and  shear  stress  distributions  change  considerably 
when  friction  forces  are  included.  The  addition  of  friction  effects 
appears  to  be  the  key  in  understanding  the  stress  distribution  in  con¬ 
tact  interfaces  between  ceramic  and  metal  components,  especially 
blade/disk  interfaces.  Therefore,  classical  solutions  for  bodies  in 
contact  with  both  normal  and  tangential  loads  were  examined.  Figure  7 
shows  the  tangential  stresses  for  various  ratios  between  the  normal 
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Figure  1.  original  Finite-Element  Model 
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Figure  2.  Refined  Finite~Eleaent  Model. 
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Figure  3.  Final  Zoom  Finite-Element  Model. 
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Figure  4.  Relative  Size  of  Pinite-Bleoent  Model. 
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Figure  6.  Comparison  of  Finite  Element  to  Closed-Form  Solution 


Figure  7.  Contact  Tangential  Stresses  for  Different  "f"  Value 
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and  tangential  loads.  In  the  absence  of  tangential  load,  the  stresses 
are  compressive  or  non-existent  along  the  surface.  With  the  addition 
of  tangential  load,  a  tensile  stress  is  generated  at  the  edge  of  the 
contact  zone.  The  tensile  load  increases  as  the  tangential  load  is 
increased  as  long  as  there  is  no  relative  motion  between  the  two 
bodies. 

The  peak  tensile  stress  is  very  localized  at  the  surface  as  well 
as  at  the  edge  of  the  contact  zone.  Figure  8  illustrates  this  stress 
peak  at  the  surface.  The  peak  tensile  stress  is  considerably  lower  in 
magnitude,  0.001  centimeter  below  the  surface. 

Zoom  modeling  is  required  to  determine  the  localized  peak  tensile 
stresses  for  contact  interface  problems  using  finite-element  computer 
programs.  Figure  9  compares  a  finite-element  solution  utilizing  zoom 
modeling  to  a  closed-form  solution  at  the  interface  for  a  cylinder  in 
contact  with  a  flat  plate  for  a  normal/tangential  load  ratio  equal  to 
3.  Seven  levels  of  zoom  modeling  were  incorporated  into  the  finite- 
element  model  to  determine  the  peak  tensile  stress.  The  zoom  finite- 
element  mesh  for  determining  the  peak  tensile  stress  is  illustrated  in 
Figure  10.  Figures  9  and  11  show  the  localized  nature  of  the  peak 
tangential  tensile  stress. 

2.2  Empirical  Approach  and  Accomplishments 

In  addition  to  the  analytical  objectives  of  this  program,  speci¬ 
men  testing  and  data  correlation  are  required  to  determine  the  failure 
mechanism  and  to  formulate  a  failure  theory  for  the  complex  state  of 
stress  in  the  contact  interface  zone.  This  verification  testing  was 
accomplished  by  simultaneously  applying  an  axial  load  and  a  compres¬ 
sive  normal  contact  load  to  the  ceramic  test  specimens,  Figure  12* 
The  test  fixture.  Part  No.  STE-255277,  is  Illustrated  in  Figure  13. 
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Figure  8.  Contact  Tangential  Stresses  Below  the  Surface. 
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-  Axial  Load 

F2  **  Compressive  Normal  Load 


Figure  12.  Ceramic  Teat  Fixture  Concept 


Figure  13.  Ceramic  Specimen  Test  Fixture. 
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Hot-pressed  silicon  nitride  rectangular  specimens  (3  in.  x  1/4 
in,  X  1/8  in.)  were  tested  in  contact  with  hardened  4340  steel  cylin¬ 
ders.  One  of  the  cylinders  was  2.0  inches  in  diameter  and  the  other 
cylinder  was  3.0  inches  in  diameter.  The  two  cylinders  were  not 
identical  to  force  the  fractures  to  initiate  on  only  one  surface  of 
the  specimen.  The  stresses  caused  by  the  contact  loads  were  approxi¬ 
mately  25  percent  higher  on  the  surface  in  contact  with  the  smaller 
cylinder. 

The  first  seven  samples  were  instrumented  with  strain  gauges  on 
both  sides  of  the  specimen  to  evaluate  the  amount  of  bending  present 
during  testing.  All  the  fractures  from  this  first  set  of  samples 
initiated  on  the  surface  of  the  smaller  diameter  cylinder.  Based  on 
these  results,  it  was  deemed  necessary  to  only  instrument  the  side  of 
the  specimen  in  contact  with  the  smaller  diameter  cylinder. 

Testing  was  conducted  by  first  placing  a  compressive  contact  load 
on  the  specimen.  The  axial  load  was  then  increased  until  slipping  or 
fracture  occurred.  The  compressive  contact  load,  tensile  axial  load, 
and  the  axial  strain  on  the  specimen  side  near  the  smaller  diameter 
cylinder  were  measured.  Testing  was  also  conducted  with  HS2S  and 
platinum  compliant  material  between  the  contacting  bodies. 

The  coefficient  of  friction  data  exhibited  good  correlation  with 
relatively  small  variations  between  test  specimens.  The  average  coef¬ 
ficient  of  friction  between  the  hardened  4340  steel  and  the  hot- 
pressed  silicon  nitride  was  0.143.  The  measured  data  varied  from 
0.113  to  0.170.  Figures  14  and  15  are  two-  and  three-parameter 
Meibull  plots  of  the  coefficient  of  friction  data,  which  illustrates 
its  well-behaved  nature.  These  data  agree  well  with  the  published 
results  for  H50  tool  steel  on  hot-pressed  silicon  nitride,  Reference 
(3). 
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Fracture  data  were  accumulated  from  1500  to  4000  pounds  compres¬ 
sive  load.  In  order  to  correlate  and  normalize  all  these  data,  the 
stresses  due  to  the  contact  load  were  superimposed  on  the  stresses 
caused  by  the  axial  tensile  load.  A  classical  solution.  Reference 
(1) ,  was  used  to  calculate  stresses  in  each  specimen  due  to  the  con¬ 
tact  loading.  These  calculated  stresses  were  superimposed  on  the 
axial  tensile  stresses  to  determine  the  maximum  tensile  and  maximum 
principal  stresses  at  failure  or  when  the  specimen  slipped. 

The  resulting  stresses  for  both  the  fractures  and  slips  are 
illustrated  through  the  use  of  two-  and  three-parameter  Weibull  plots. 
Figures  16  and  17.  The  slip  data  are  included  in  the  statistical 
r^nalysis  as  suspended  items  and  the  fracture  data  as  failures.  Based 
on  the  two-parameter  Weibull  plot,  the  characteristic  strength  is  94.5 
KSI  and  the  Weibull  modulus  is  4.2.  Figure  18  compares  the  contact- 
f cacture-strength  results  to  four-point  and  three-point  flexure  data 
acquired  from  testinq  of  specimens  from  the  sane  billet  of  material. 
Under  ideal  test  conditions,  the  resulting  characteristic  strength 
would  be  expected  to  be  better  than  four-point  flexure  results  and  to 
be  similar  to  three-point  flexure  results  because  the  amount  of 
stressed  surface  area  during  these  tests  is  closer  in  size  to  the 
stressed  surface  area  for  three-point  flexure  tests.  The  scatter  in 
the  data  can  probably  be  attributed  to  specimen-to-specimen  variations 
in  surface  finish  in  the  interface  zone,  alignment  of  the  specimen  In 
the  test  fixture,  and  alignment  of  the  contacting  cylinders.  As  men¬ 
tioned  previously,  the  peak  strese^s  are  extremely  localized,  even 
more  so  than  three-point  flexure  stresses.  These  localized  stresses 
could  possibly  be  affected  by  a  different  distribution  of  flaws  than 
flexure  stresses. 

Another  possibility  for  introducing  errors  into  the  fracture 
strength  results  are  the  stress  calculations  accounting  for  the  con¬ 
tact  loading.  These  two-dimensional  stress  calculations  assume  con¬ 
stant  stress  along  the  entice  width  of  the  test  specimen.  In  actu¬ 
ality,  these  stresses  are  not  constant.  Even  ignoring  misalignment  of 
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the  contacting  cylinders,  the  state  of  stress  at  the  edge  of  the 
specimen  is  plane  stress  and  at  the  middle  of  the  specimen  is  plane 
strain. 

The  test  results  are  very  reasonable  and  verify  the  maximum  prin¬ 
cipal  stress  theory  for  failure  prediction  in  brittle  materials.  All 
the  fractures  which  have  been  examined  to  date  indicate  that  the  frac¬ 
tures  originated  on  the  surface  at  the  edge  of  the  contact  zone,  which 
is  the  location  of  maximum  principal  stress. 

Testing  has  also  been  completed  using  HS25  and  platinum  compliant 
materials  between  the  contacting  surfaces. 

Coefficient  of  friction  data  for  these  two  compliant  layers  on 
hot-pressed  silicon  nitride  is  illustrated  in  Figure  19.  The  specimen 
could  tolerate  higher  contact  loads  when  the  compliant  layers  were 
utilized.  The  higher  loads  were  especially  evident  when  the  platinum 
layers  were  tested.  However,  the  platinum  was  so  ductile  during  test¬ 
ing  that  fracture  strengths  were  impossible  to  measure.  The  platinum 
simply  extruded  until  it  was  too  thin  to  be  an  effective  compliant 
medium.  Consequently,  in  order  to  evaluate  the  life  capabilities  of 
"soft**  compliant  media  similar  to  platinum,  it  is  necessary  to  perform 
cyclic  specimen  tests. 

Typical  fracture  surfaces  are  shown  in  Figure  20.  Some  of  the 
fracture  origins  tended  to  be  diffused  and  were  not  easily  identified. 
Two  such  cases  ate  shown  in  Figures  20(a)  and  20(b).  In  other  cases, 
Figure  20(c),  the  origin  could  be  identified  but  was  not  as  distinct 
as  is  typical  of  flexural  test  specimens  of  hot-pressed  i  Ilicon 
nitride.  This  characteristic  is  probably  related  to  the  localized 
nature  of  the  contact  stresses,  and  the  relatively  small  amount  of 
stored  elastic  energy  at  fracture. 
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(a)  SPECIMEN  NO.  18 


(22X 


(12X) 

(b)  SPECIMEN  NO.  14 


Figure  20.  S£M  Photographs  of  Fracture  Surfaces 


(12X)  (25X) 

(d)  SPECIMEN  NO.  8  (ARROW  INDICATES HERTZIAN  CONE  CRACK) 

Figure  20.  (Contd.)  SEM  Photographs  of  Fracture  Surfaces. 
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In  several  cases,  there  were  indications  that  fracture  initiated 
from  a  classical  Hertzian  cone  crack.  An  example  of  fracture  from  a 
cone  crack  is  shown  in  Figure  20  (d) ,  This  type  of  crack  is  due  to  a 
contact  stress  localized  at  a  point  —  possibly  due  to  a  speck  of  dirt 
between  the  cylindrical  contact  and  the  specimen. 

Another  feature  of  the  fracture  surface  which  differed  from 
specimen  to  specimen  was  the  topography.  Some  specimens.  Figures 
20(c)  and  (d) ,  had  relatively  smooth  fracture  surfaces  while  others. 
Figures  20(a)  and  (b)  ,  exhibited  regions  of  considerable  curvature. 
The  exact  reason  for  this  difference  is  not  known,  but  should  be 
related  to  the  contact  stress  distribution  and  the  location  of  frac¬ 
ture  initiation. 
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3.0  CONCLUSIONS 

The  following  conclusions  have  been  reached  as  a  result  of  this 
contact- stress  program. 

(a)  The  finite-element  approach  using  zoom  modeling  can  be  used 
to  predict  the  state  of  stress  in  ceramic-to-metal  contact 
interface  zones. 

(b)  The  probability  of  success  for  a  ceramic-to-metal  interface 
can  be  predicted  using  maximum  principal  stresses. 

(c)  Compliant  media  can  increase  the  allowable  contact  loads 
which  a  ceramic  can  tolerate. 

(d)  More  understanding  of  the  necessary  physical  properties  of  a 
compliant  media  is  required  for  a  long-life  application. 
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